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Abstract. The FIRBACK (Far InfraRed BACKground) 
survey represents the deepest extensive 170 /im images ob- 
tained by the ISO satellite. The sources detected com- 
prise about 10% of the Cosmic IR Background (CIB) seen 
by COBE, and, importantly, were observed at a wave- 
length near the peak of the CIB. Detailed follow-up of 
these sources should help pin down the redshifts, Hubble 
types, and other properties of the galaxies which consti- 
tute this background. We have used the Submillimetre 
Common-User Bolometer Array (SCUBA) instrument on 
the James Clerk Maxwell Telescope (JCMT) to search for 
sub- mm emission from a sample of 10 galaxies which are 
bright at 170 fim and for which we had accurate radio po- 
sitions. Statistically we detect this sample at both 450 fim 
(~ 4(t) and 850 /im (~ 7ct); individual objects are convinc- 
ingly detected in four cases at 850 /xm, and one case at 
450 /zm. Fits to simple spectral energy distributions sug- 
gest a range of low to moderate redshifts, perhaps z = 0- 
1.5 for this FIRBACK sub-sample. 
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1. Introduction 

In the far-IR/sub-mm waveband there are currently two 
pressing (and related) cosmological mysteries: 

1. what makes up the background detected by COBE 
(Puget et al. 1996, Fixsen et al. 1998, Hauser et 
al. 1998, Lagache et al. 1999)? 
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2. how do sources detected in the sub- mm (see 
e.g. Sanders 1999 and references therein) relate to 
galaxies at other wavelengths? 

The study of both issues provides important diagnos- 
tics of galaxy formation models, and the answers should 
help illuminate the dark ages of how the first objects 
formed and subsequently evolved into present day galax- 
ies. The first question is the main motivation behind the 
FIRBACK project (Clements et al. 1998, Lagache 1998, 
Reach et al. 1998, Puget et al. 1999, Dole et al. 1999), 
which made deep 170 /im images of the sky with ISO, 
to resolve the CIB into sources. The second question has 
been the main driving force behind the search for dis- 
tant sources by several teams using SCUBA, as well as 
the large amount of follow-up work and comparison with 
other wavelengths (e.g. Small, Ivison & Blain 1997, Barger 
et al. 1998, Blain, Ivison & Small 1998, Holland et al. 1998, 
Hughes et al. 1998, Small et al. 1998, Barger et al. 1999, 
Blain et al. 1999a, Chapman et al. 1999, Cowie et al. 1999, 
Bales et al. 1999, Lilly et al. 1999). These observational 
campaigns are untangling the problem of how obscuration 
by dust skews our optical view of the early Universe, un- 
veiling the 'dark-side' of galaxy formation out to distant 
redshifts, and helping provide unbiased estimates of the 
global star formation history of the Universe. 

In this letter we present the first results of a study 
at the interface between these two puzzles. Rather than 
trying to find new sources with SCUBA, we have carried 
out 450 fiTa and 850 /im photometry on sources already 
detected at 170 /im with ISOPHOT. 

2. The FIRBACK Survey 

The lifetime of the Infrared Space Observatory {ISO; 
Kessler et al. 1996) is long over, and we now have 
in hand the best information we will have from long- 
IR wavelengths until the launch of SIRTF. For SCUBA 
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observations, the smallest extrapolations come at the 
longest wavelengths attainable with the ISOPHOT in- 
strument (Lemke et al. 1996). ISOPHOT is an imaging 
photo-polarimeter with 92 arcsec pixels and 1.6 arcmin 
FWHM at 170 ^m; reasonably high signal-to- noise sources 
in 'dithered' images can be located with an accuracy 
of ~40 arcsec. Following the discovery of the GIB, ISO 
was in a unique position to investigate the sources which 
comprise this background. Consequently deep ISOPHOT 
images were obtained at 170 /xm of three separate '--^l 
square degree regions of the sky, selected for their low 
foreground emission - the FIRBACK Survey. One region, 
the 'Marano' fields, is only accessible from the southern 
hemisphere, while the other two coincide with the 'Nl' 
and 'N2' fields of the ELAIS project (Oliver 1997), which 
used ISOCAM at 7 ^im and 15 //m and ISOPHOT at 90 //m 
(although only a fraction of the FIRBACK sources were 
detected at these other wavelengths). These northern re- 
gions have also been mapped with the VLA (Ciliegi et 
al. 1999). In addition an area towards the Lockman Hole 
is being studied by another group (Kawara et al. 1998). 

Here we have concentrated on a sample of objects 
from the roughly two square degree 'Nl' field, centred at 
161^ ll''^, +54° 25™. For this first attempt at JCMT follow- 
up we chose objects with secure, unconfused radio identi- 
fications, relatively strong 170 /um emission, and addition- 
ally high 170/zm:21cm flux density ratios (see Lagache et 
al., in preparation). This last criterion was aimed at bi- 
assing the sample away from the lowest redshift galaxies, 
and hence towards those which might have the highest 
850/im:170/im flux density ratios. We expect that with 
the addition of more follow-up observations it should be 
possible to select future sub-samples with a higher likeli- 
hood of being strong SCUBA emitters. 

3. JCMT Observations 

The data were obtained on the nights of 18-23 March 
1999 using the SCUBA instrument (Holland et al. 1999) 
on the JCMT. The short- and long- wavelength arrays were 
used simultaneously, at 450 /xm and 850 /.tm, respectively. 
We used 'photometry' mode, chopping at the standard 
7.8125 Hz, and also nodding every second by about 45 
arcsec in coordinates fixed to the array (i.e. there was 
no sky rotation). This means that each measurement is a 
double-difference between the central bolometer and po- 
sitions 45 arcsec each side, corresponding approximately 
to positions of other bolometers on both the long- and 
short-wavelength arrays. The data were analyzed using 
the SURF package (Jenness & Lightfoot 1998). The raw 
data for the two arrays were flat-fielded, corrected for ex- 
tinction, had bad bolometers removed, and had the aver- 
age sky removed at each time interval. The information 
from the off-beams was then added, assuming that one 
long wavelength bolometer had an efficiency of exactly 
0.5, with the long wavelength bolometer on the other side. 
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Table 1. FIRBACK identification (see Lagache et al., 
in preparation), together with 170 /^m flux density (with 
estimated systematic uncertainty of perhaps 30%). The 
450 fim and 850 fiia measurements are from our new 
SCUBA observations. The final column lists the SNR for 
the > 3ct detections, and the 95% Bayesian upper limits 
(in mJy) for the others. 



as well as the two short wavelength off-beam bolometers, 
having slightly lower values (see Borys et al. 1999, Chap- 
man et al. 1999, for more details). Adding the weighted 
off-beam signal always decreased the noise, and generally 
increased the signal to noise ratio (SNR), but we carried 
out the same procedure even when it slightly lowered the 
SNR. 

Calibration was performed a few times per night using 
planets and other strong sub-mm sources, and the values 
we used were similar to the standard gains. The standard 
deviation of the calibrations was 10% at 850 /^m and 12% 
at 450 fim; these should be a reasonable estimate of the 
uncertainty in the calibration. At the low SNRs at which 
we are working, the calibration uncertainty is not a major 
contributor to the total uncertainty, and has essentially 
no effect on the SNR itself. 



4. Individual Objects 

At 850 ^m we detected one source at > 5a and a further 
three at > 3a. While we would not claim detection of the 
other sources, they generally have positive flux density 
(see next section) , and certainly there are good upper lim- 
its in each case. This last remark applies to the 450 /xm 
data, where there is only one detections above 3ct. 

Bayesian 95% upper limits can be obtained for all our 
non-detections, by integrating a Gaussian probability, ne- 
glecting the unphysical negative flux density region. The 
850 fim upper limits for our six non-detections are given in 
the last column of Tabic ^ At 450 /im the limits are gen- 
erally less constraining for reasonable SEDs, being around 
30 mJy in the best cases. 

Using a combination of 170 /im, 450 /im and 850 /xm 
data, we can place some constraints on the spectral en- 
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ergy distributions (SEDs) of these FIRBACK galaxies. 
Assuming a grey-body spectrum, we can obtain a limit on 
sonic combination of luminosity, temperature, spectral in- 
dex and redshift. Here we choose to normalize the luminos- 
ity at 170 fim, and then use the SCUBA data to constrain 
the redshift. To do this we assume standard values for the 
dust temperature, Td = 40 K (typical for sub-mm selected 
galaxies, e.g. Blain et al. 1999b, Dunne et al. 1999), and 
spectral index of the dust cmissivity, (3 — 1.5. Because we 
do not know the absolute luminosity of any source, our 
results are degenerate in the ratio Td/(1 -I- z), and so we 
are unable to tell apart cooler objects at lower redshift 
from hotter objects at higher redshift. 

Assuming a uniform prior distribution of redshifts, we 
can obtain a Bayesian 95% confidence range on the red- 
shift implied for each source by our 450 /im and 850 /im 
data, where it is understood that the redshifts can be 
scaled using a different dust temperature such that (1 + 
z)/Td is held constant. We find that the objects with the 
highest 850 /im flux densities have the highest implied red- 
shifts: e.g. 0.66<z<1.23 for Nl-038. Those with lower 
flux densities at 850 /im generally only yield upper limits 
to the redshift: e.g. z < 0.55 for Nl-015. There are no ob- 
jects for which we would infer z > 2. However, higher red- 
shifts could be accommodated by adopting a higher dust 
temperature or a higher value of /? (and lower redshifts by 
lowering these parameters) . 

The object in our sub-sample which is brightest at 
170 /im, Nl-008, is hard to flt with any reasonable SED; 
it has approximately half the 850 /xm and 450 /^m flux 
density expected from even a z = source with the same 
170 /im flux density. This suggests that there might be 
more than one source contributing to the ISOPHOT flux, 
which is not particularly unlikely, since the FIRBACK 
Survey is operating near the confusion limit. On the other 
hand, optical images show little in the error circle except 
for a very bright obvious spiral galaxy, and so from that 
point of view this is not a case where we expect more than 
one source in the ISOPHOT beam. Another possibility is 
that our JCMT beam did not include all the flux, since 
optically this object appears quite extended. However, we 
would expect the sub-mm emission to be more concen- 
trated than the optical, and hence we are likely to have 
included most of the dust emission within the beam. This 
object will be discussed more extensively in Lagache et 
al. (in preparation). 

Relatively poor values of x^ are also found for Nl-034 
and Nl-063. These result either from low SCUBA flux 
densities relative to ISOPHOT, or somewhat high SCUBA 
450 /im flux densities compared with 850 /im. Higher SNR 
data, or data at additional wavelengths are required to 
determine whether these flux densities come from multiple 
objects, complex SEDs, or simply the low SNRs of the 
current data for these objects. 

Another way of estimating the redshift relies on the 
radio/far-IR correlation (e.g. Helou, Soifer & Rowan- 
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Fig. 1. Average spectral energy distribution for our FIR- 
BACK sub-sample. The 170 /im point is the average of 
our sub-sample of 10 FIRBACK sources, with the error- 
bar here being the standard error from the scatter among 
them. The 450 /im and 850 /im points are from our new 
SCUBA observations. For the sake of visual clarity we 
have multiplied the y-axis by x^. The curves show emis- 
sion from modified blackbodies, normalized to the 170 /im 
flux density, with T^ = 40K and /3=1.5, and for z = 0.0 
(solid line) up to z = 1.0 in steps of 0.2. Note that the 
shapes of these curves are degenerate in the combination 
{l + z)/Td. 

Robinson 1987, CariUi & Yun 1999, Barger, Cowie & 
Richards 1999, Smafl et al. 1999). Using the 20 cm VLA 
data from Ciliegi et al. (1999) and the explicit correla- 
tion (using (3 = 1.5, Td = 40 K) from Carilh & Yun (1999), 
we find z = 0.6, 1.2, 1.1 and 1.4 for Nl-008, Nl-038, Nl- 
061 and Nl-063, respectively (our four 850 /zm detections). 
Except for Nl-008 these are in broad agreement with the 
values obtained from the sub-mm and far-IR data alone. 

5. Statistical Results 

If we combine all the data together, then we can ob- 
tain a much more precise picture of the average galaxy 
in our FIRBACK sub-sample. The average flux densities 
are given in Table 1, and represent a 4.1(7 detection at 
450 /im and a 7.3cr detection at 850 /im. The average SED 
is shown in Figure 1. It is clear that for Td = 40 K the best 
fitting redshift is around 0.3, with values as low as z = or 
higher than z ~ 0.6 providing relatively poor fits. However, 
we can certainly accommodate z ~ 1 galaxies for higher 
temperatures, and z ^0 galaxies for lower temperatures. 
Since higher dust temperature is seen in star-bursting 
galaxies, compared with more typical star-forming galax- 
ies, then it is possible that a fraction of these sources are 
much more luminous and at higher redshift. However, it 
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would take unrealistically high temperatures, T^ ^ 100 K, 
to push some of the objects up to say z ~ 3 (although the 
possibility of gravitational lensing could complicate this). 
Models of galaxy populations (e.g. Guiderdoni et 
al. 1998), designed to fit number counts at a range of 
wavelengths, as well as the GIB, predict that the average 
FIRBAGK galaxy is indeed at z ~ 1. The other possibility 
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is thati wti have uiicovyred a new pupulttLiuii of rultiLivdy 
nearby sLai'-fui'iiiiiig galaxies, which do iiuL show up clearly 
in surveys at other wavelengths. The degeneracy between 
(1 + z ) and 1/Td makes it impossible to decide between 
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these p pssibilitics using the far-IR and sub-mm data alone. 
Although it seems unlikely that all our objects can be at 
low redshift, this could obviously be resolved by obtain- 
ing optical redshifts. Once redshifts have been obtained, 
then t he oub mm data will bo invaluable in moaouring the 
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proper tioo of the duot in the variouo gala^cy typoo: far IR 
luminosities, dust temperatures and emissivities. 

6. Conclusions 

We have carried out the first SGUBA follow-up of FIR- 
BAGK sources. We found that they are generally de- 
tectable in the sub-mm; those with somewhat higher 
850 /xm flux density may be at z~l, while those which 
are fainter in the sub-mm may be more normal galax- 
ies at z r^O. Models of evolving galaxy populations which 
provide a good fit to the 170 /im counts, as well as counts 
at other wavelengths (e.g. Guiderdoni et al. 1998) pre- 
dict that the median redshift of the FIRBAGK galaxies is 
around 1 (Puget et al. 1999). Our results are consistent 
with this, provided that the average galaxy in our sub- 
sample is a distant star-bursting galaxy with fairly hot 
dust temperature. The other possibility is that some could 
be from an otherwise unknown population of low redshift 
star-forming galaxies with relatively low dust tempera- 
ture. Further observations at sub-mm and other wave- 
lengths should decide this issue, and reveal the detailed 
properties of the galaxies which comprise the GIB. 
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